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It  is  concluded  that  provided  the  electrical  properties  of  compacted  lead 
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Abstract 

Lead  telluride  type  semiconductors  are  used  in  the  fabrication  of 
thermoelectric  modules  ^currently  employed  in  a  number  of  US  military 
applications. -* This  report  covers  the  programme,  of  research  undertaken  in 
the  Department  of  Physics,  Electronics  and  Electrical  Engineering  at  UWIST 
Cardiff  'to  produce  materials  based  upon  lead  telluride  with  improved 
"figures  of  merit" '  and  hence  greater  thermoelectric  conversion 
efficiency. 

One  way  of  improving  the  figure  of  merit  is  by  reducing  the  lattice 
thermal  conductivity  of  the  material.  This  can  be  achieved  by  increasing 
phonon^grain  boundary  scattering.  A  realistic  theoretical  model  has  been 
developed  for  lead  telluride  and  used  to  investigate  the  lattice  thermal 
conductivity  as  a  function  of  grain  ,  size  and  level  of  doping.  In f 
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Phonon  grain  boundary  scattering  is  enhanced  in  semiconductor  alloys 
because  of  the  presence  of  disorder  scattering  and  the  theoretical  model 
was  extended  to  take  this  factor  into  account.  PbSnTe  and  PbGeTe  were 
identified  as  alloys  whose  lattice  thermal  conductivity  could  be 
significantly  decreased  by  a  reduction  in  grain  size  and  in  optimaly  doped 
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equivalent  single  crystal  material  was  estimated  to  be  11  and  14  percent 
respectively. ^ 

It  is  concluded  that  provided  the  electrical  properties  of  compacted 
lead  telluride  material  can  be  maintained  close  to  that  of  single  crystal; 
the  thermoelectric  figure  of  merit  can  be  substantially  improved  by 
employing  a  fine  grain  size.  In  PbGeTe  with  a  grain  size  of  1/im,  the 
improvement  in  the  thermoelectric  figure  of  merit  is  estimated  to  be  about 
10  percent. 
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Table  1. 


Important  material  parameters  for  lead  tellunde. 


Table  2. 

Table  3. 


Xe/VL  f°r  ■lead  telluride  corresponding  to  $opt.  (non¬ 
parabolic  bands,  acoustic  phonon  scattering. 

xe/vL  for  lead  telluride  corresponding  to  {opt.  (optical 
phonon,  acoustic  phonon  and  combined). 


II  GENERAL  INTRODUCTION 


Multifuel  thermoelectric  generators  employing  modules  based  upon  lead 
telluride  technology  are  used  in  a  number  of  military  applications1 . 

In  addition  to  their  multifuel  capability  these  generators  are  difficult 
to  detect  by  visual,  aural,  or  thermal  infrared.2  Clearly  any  improvement 
in  the  thermoelectric  conversion  efficiency  of  the  generator  would  result 
in  a  saving  in  fuel  -  an  important  consideration  when  the  device  is  used 
in  a  tactical,  situation. 

The  factors  which  determine  the  conversion  efficiency  (n)  of  a 
thermoelectric  generator  have  been  exhaustively  dealt  with  in  the 
literature.3"11  It  is  suffice  to  state  that  once  the  operating 
temperature  and  the  temperature  difference  over  which  the  generator 
operates  has  been  decided,  n  depends  solely  upon  the  so  called  figure  of 
merit  Z  of  the  thermocouple  material,  where  Z  =  <x2o/x,  a  is  the  Seebeck 
coefficient,  o  the  electrical  conductivity.  The  thermal  conductivity  x  is 
the  sum  of  xL,  a  contribution  due  to  the  lattice  and  xe  an  electron  (or 
hole)  contribution. 

The  parameters  which  occur  in  the  expression  for  the  figure  of  merit 
are  functions  of  the  carrier  concentration  (usually  expressed  in  terms  of 
the  reduced  Fermi  energy  ^ )  and  in  the  established  high  temperature 
thermoelectric  material  silicon-germanium  alloy,  Z  is  optimized  at  around 
102  5  -I02  6m"3  12 .  However,  at  these  high  carrier  concentrations  xL 
still  accounts  for  about  75  percent  of  x  and  it  is  anticipated  that  xL 
will  contribute  a  similar  proportion  to  the  thermal  conductivity  of  lead 
telluride.  Consequently  a  reduction  in  xL  would  result  in  a  corresponding 
increase  in  the  figure  of  merit  of  lead  telluride  and  hence  in  the 
materials  thermoelectric  conversion  efficiency. 

It  has  been  reported  that  phonon-grain  boundary  scattering  has  a 
significant  effect  on  reducing  the  lattice  thermal  conductivity  of  silicon 
germanium  alloys.13-15  This  phenomenon  does  not  appear  to  be  accompanied 
by  a  deterioration  in  the  other  parameters  which  occur  in  the  figure  of 
merit  although  the  long  term  behaviour  of  small  grain  size  material 
when  subjected  to  high  temperatures  has  not  been  reported.  Consequently 
small  grain  size  silicon  germanium  alloys  exhibit  a  higher  figure  of  merit 
than  comparable  "single  crystal"  or  large  grain  size  material.  Grain 
boundary  scattering  is  particularly  favoured  in  silicon-germanium  alloy 
because  the  large  difference  in  atomic  masses  of  the  constituent  atoms 
give  rise  to  substantial  alloy  disorder  scattering.  However,  this 


phenomenon  will  also  be  present  in  other  thermoelectric  alloys  and  in 
alloy  compounds  such  as  those  based  upon  lead  telluride. 

The  United  States  Army  is  currently  evaluating  the  performance  of  a 
number  of  improved  thermoelectric  materials  for  possible  use  in  future 
thermoelectric  generators.  Evidently  information  on  the  reduction  in 
thermal  conductivity  and  hence  potential  improvement  in  thermoelectric 
performance,  which  accompanies  the  use  of  small  grain  size  lead  telluride 
type  material  and  on  the  conditions  which  optimise  the  thermoelectric 
figure  of  merit  of  these  materials  is  relevent  to  their  programme  of 
material  evaluation. 

III  OBJECTIVES 

The  objectives  of  the  programme  of  research  described  in  this  report 
were : - 

1(a)  Initially  to  develop  a  semiquantitative  theoretical  model  of  lead 
telluride.  Use  the  model  to  estimate  the  relative  reduction  in  the 
lattice  thermal  conductivity  compared  to  single  crystal  material  with 
decrease  in  grain  size  and  hence  improvement  in  the  thermoelectric 
figure  of  merit. 

(b)  Substantiate  the  theoretical  prediction  by  preparing  <5jxm  grain  size 
compacts  of  lead  telluride  and  measuring  the  relevant  transport 
properties. 

2  Depending  upon  the  success  or  otherwise  in  achieving  these  initial 
objectives  the  second  phase  of  the  work  was  to:- 

(a)  Develop  a  realistic  theoretical  model  for  the  absolute  magnitude  of  the 
thermal  conductivity  of  lead  telluride.  Use  the  model  to  identify 
alloys  based  on  lead  telluride  whose  lattice  thermal  conductivity 
offer  the  best  potential  for  reduction  as  a  result  of  phonon-grain 
boundary  scattering.  Determine  the  conditions  for  optimising  the 
thermoelectric  figure  of  merit  of  these  alloys. 

(b)  Attempt  to  comminute  very  small  grain  size  (<0.5>im)  lead  telluride 
material . 

For  convenience  in  presentation  the  work  reported  is  divided  into  two 
stages  1.  Theoretical  and  2.  Experimental 

IV  DEVELOPMENT  OF  A  THEORETICAL  MODEL  FOR  LEAD  TELLURIDE 
1.  Introduction 

The  initial  objective  to  estimate  the  relative  reduction  in  lattice 
thermal  conductivity,  compared  to  single  crystal,  with  decrease  in  grain 
size  was  fairly  straightforward.  A  two  band  model  with  parabolic 
multivalued  energy  band  structure  was  considered.  Acoustic  phonon 


scattering  was  taken  to  be  the  dominent  scattering  mechanism,  intervalley 
scattering  is  neglected  and  no  distinction  made  between  the  conductivity 
effective  mass  and  the  density-of -states  effective  mass.  Although  the 
model  was  later  improved,  preliminary  calculations  indicated  that  the 
relative  changes  in  thermal  conductivity  due  to  phonon-grain  boundary 
scattering,  and  consequent  relative  changes  in  the  thermoelectric  figure 
of  merit  are  fairly  insensitive  to  the  inclusion  of  these  refinements. 

However,  the  development  of  a  realistic  model  for  lead  telluride  for 
use  in  identifying  potentially  good  alloys  based  on  lead  telluride 
presented  considerable  difficulties.  The  theoretical  model  must  give  good 
agreement  with  experimental  data  cited  in  the  literature  and  this  involved 
obtaining  the  separate  electronic  (xe)  and  lattice  (xL)  contributions  to 
the  thermal  conductivity.  A  prerequisite  for  this  analysis  is  a  knowledge 
of  the  Lorenz  number.  A  literature  survey  revealed  that  there  is  very 
little  European  or  United  States  published  information  on  the  value  and 
behaviour  of  the  Lorenz  number  for  lead  telluride.  The  realistic  model 
developed  considers  both  acoustic  and  optical  phonon  scattering.  It  takes 
into  account  the  non  parabolic  nature  and  multivalued  structure  of  the 
energy  bands,  including  intervalley  scattering.  Minority  carrier  effects 
will  also  be  significant  over  part  of  the  temperature  range  of  operation 
of  the  materials  and  should  be  included.  The  principal  elements  of  the 
theoretical  model  are  discussed  in  the  following  section. 

2.  Lattice  thermal  conductivity 


At  temperatures  above  room  temperature,  a  simple  formulation  for  the 
lattice  thermal  conductivity  can  be  used  which  includes  the  various  phonon 
scattering  mechanisms16'17.  The  lattice  thermal  conductivity  is  expressed 
in  terms  of  three  parameters  A,  B  and  C  which  relate  to  phonon  scattering 
by  alloy  disorder,  free  carriers  and  grain  boundaries  respectively. 
Defining  xQ  as  the  thermal  conductivity  of  a  large,  perfect  (notional) 
crystal,  the  lattice  thermal  conductivity  is  given  by:- 
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where  L^AjBfC)  = 


Ax4+x2+8x+C 


here  x  =  fiu/kBT  and  ko  =  Tu/TN  *s  the  rat*0  of  the  phonon  relaxation  times 
for  Umklapp  and  Normal  processes,  kB  is  Boltzmann’s  constant  and  T  is  the 
temperature.  There  is  no  advantage  in  this  present  study  to  undertake 
detailed  calculations  of  the  contributions  to  the  thermal  conductivity 
from  separate  polarizations  and  an  average  acoustic  branch  is  assumed. 
The  various  phonon  relaxation  times  are  as  described  in  the 
literature;14'16'17  the  phonon-electron  scattering  relaxation  time  is 
given  by 
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The  various  parameters  have  the  same  meaning  as  in  the 
literature14'1617.  It  is  usual  to  express  C  in  terms  of  a  parameter  D, 
which  is  inversely  proportional  to  the  grain  size  L  and  they  are  related 


by  D=CT  where  T  is  the  temperature.  A=0  corresponds  to  unalloyed  material 
with  no  disorder  present,  B=0  corresponds  to  undoped  material  and  C=0  to 
single  crystal  material.  In  general  x(A,B,C=0)  =  ^singie  represents  the 
lattice  thermal  conductivity  of  a  doped  single  crystal  alloy  while 
x  (A,B,C)  =  ^sintered  represents  the  lattice  thermal  conductivity  of  a 
compacted  (sintered)  alloy. 

3.  Electronic  contribution  to  the  thermal  conductivit 


(1)  Introduction 

Lead  telluride  and  its  alloys  possess  a  narrow  energy  band  gap  (Eg). 
These  materials  have  small  electron  effective  masses  and  the  density  of 
states  which  varies  as  m *3A  is  also  relatively  small.  Consequently  a 
relatively  small  number  of  carriers  will  be  sufficient  to  fill  the  bands 
up  to  high  energy  levels.  The  effective  mass  of  electrons  near  the  top  of 
the  conduction  band  become  energy  dependent18*  In  the  model  adopted  a  two 
band  conduction  band  is  considered  with  the  following  main  features;  the 
band  extrema  for  the  conduction  and  valence  bands  are  assumed  to  be 
located  at  the  same  k  value.  The  energy  separation  from  the  other  bands 
at  this  k  value  is  greater  than  the  main  energy  gap  and  the  momentum 
operator  has  non-zero  matrix  elements  between  the  states  corresponding  to 
the  extremal  points.  Non-parabolicity  and  the  effective  mass  of  both 
bands  are  defined  by  the  interaction  of  the  electron  and  hole  bands  only. 
The  energy  dispersion  law  is  of  the  form:- 


ft:kjL 


ft2k  n 


=  E  (1  +-  ) 

F 

Lg 


k,  and  k.  are  the  transverse  and  longitudinal  components  of  the  wave 

1  *  * 

vector  and  m|0  and  m||0  sstir  the  components  of  the  effective  mass  tensor 
near  the  band  extremum.  This  is  usually  referred  to  as  the  "Kane  model". 
The  energy  dependence  of  the  effective  mass  is  given  by 


*  f  2E 
n*  1+  — 
0  E 


V>\ 


•*.  - 
V-& 


v  v  V 


For  a  multivalued  structure  with  Nv  equivalent  valleys  the  density  of 
states  effective  mass  iss¬ 
u'd  =  NvV3  (»*  m^*2]  =  NvV»  nj, 

and  the  conductivity  effective  mass:- 

=  t(®  n*_i  +  ^t'1] 

The  energy  dependence  of  the  effective  mass  manifests  itself  in  the 
various  transport  coefficients;  these  are  expressed  in  terms  of  averaged 
quantities  such  as  <r>,  <r2>,  <r(E)>  etc.  In  the  case  of  non-parabolic 
energy  bands 

f  (-3f/3E)[T(E)/m*(E)]nEin 
<TnEm>  =  Jo  - 

j“  (-3f/3E) (k3/mc*n)dE 
where  k  =  (2mJ()V2  [E(l+E/Eg)]^ 

~K 

In  lead  telluride  at  room  temperature,  the  contribution  of  impurity 
scattering  to  the  carrier  mobility  is  very  small  and  will  be  negligible  at 
the  temperatures  of  operation  of  a  thermoelectric  device.  Consequently 
only  electron-phonon  scattering  mechanisms  need  to  be  considered. 

3  (ii)  Carrier  scattering  by  acoustic  phonons 

The  carrier  relaxation  time  can  be  written  as:- 

1 

T~  iMi2p(n) 

where  n  =  E/kBT  is  the  reduced  carrier  energy,  M  is  the  matrix  element  for 
the  electron-phonon  interaction  and  p(n)  is  the  density  of  states.  The 
calculation  is  simplified  if  it  is  assumed  that  the  matrix  element  depends 
upon  energy  as  for  the  case  of  a  parabolic  band.  For  scattering  by 
acoustic  phonons,  this  amounts  to  iMi2  =  constant.  The  relaxation  time  is 
then  given  by  :- 

1  1 

T - a  - - - 

ac  P(n)  [n( 1+CgH) ]”( l+23gn) 


where  6g  =  kgT/Eg 


The  appropriate  expression  for  t (n)  which  includes  the  energy  dependence 
of  M  is  given  by19'20 

nkBT(p/Nv)  f  80  n(l+B  n)] 

T(rj)-1  =  —  £2  1-  ^  * _ 

AC,  l  3(l+20gn)2 

where  Ct  is  the  elastic  constant  related  to  the  average  sound  velocity, 
(p/Nv)  is  the  density  of  state  in  a  given  valley  and  «  is  the  deformation 
potential  constant. 

The  Lorenz  factor  is  given  by:- 

£  -  (’LljAlJ- 


where  5  =  and  nL^  are  generalised. 

Fermi  integrals  and  defined  as:- 


-  —  nn[n(l+0gn)]ra(l+20gn)idn 


The  reduced  electrical  conductivity  is  given  by:- 


o’  =  kN  °Li 


v  “-z  - i - 

(®c  xl) 


With  the  carrier  concentration  n  related  to  the  reduced  Fermi  energy 
5,  by  the  expression:- 


_<X  W4  ,4 

L0 


n  = 


3n*fi 


2 « 3 


The  carrier  mobility  <iac  due  to  acoustic  phonon  scattering  is  given 


by:- 


constant 


%c 


®J0md03/z  T  3/2 


Oy'  2 

Lo 


3.{iii)  Carrier  scattering  by  polar  optical  phonons 

The  carrier  relaxation  time  for  scattering  by  polar  optical  phonons 
is  given  by20 


2I/2e2kBTmd*"2 

h2(i7kBT)ly^ 


-i 

0 


) 


l+2Bgn 

(1+V^ 


x 


[l-5„in(l+s00"1)]- 


2Bg>7  (  1+Bgfl  ) 
(l+2BgP)2 


[l+26a>+2  602-in(l+60(1  *)] 


The  Lorenz  factor  is  given  by 


where  6  =  ‘Li2/°Li2 

In  the  calculations,  screening  effects  have  not  been  taken  into  account. 
Polar  optical  scattering  although  appreciable  in  lead  telluride  materials, 
is  not  the  dominant  scattering  mechanism.  Consequently  neglecting 
screening  effects  will  not  substantially  affect  the  overall  electronic 
transport  properties. 

The  reduced  electrical  conductivity  can  be  expressed  as:- 
o'  -  NvK'T2f(0o/T)md*°Li2/(XLmc*) 


4kB2 

with  K1  =  - 

3e2n2fi(fc7— €0-1 ) 

where  e0  and  £«,  are  the  static  dielectric  constant  and  the  dielectric 
constant  at  high  frequencies  respectively.  0O  is  the  temperature  which 
corresponds  to  the  optical  phonon  frequency  at  the  zone  boundary. 

The  carrier  mobility  j/0p  due  to  optical  phonon  scattering  is  given 

by:- 


^op 


K  "  f  ( 0O  /T )  °Lf 

^  m*4°Lol 


2 


„  8ft2 (2nkB) 

where  K  =  - jj - 

3E{2nkB)  /2(*."1-*0-1) 
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3  (iv)  Multivalley  energy  band  structure  and  intervalley  scattering 

A  multi valley  energy  band  structure  together  with  intervalley 
scattering  must  be  included  in  the  theoretical  model  for  lead  telludide 
type  materials  if  quantitative  correspondence  is  sought  between  theory  and 
experimental  data.  Although  the  inclusion  of  a  multivalley  structure  and 
non-parabolicity  of  the  energy  momentum  relationship  makes  the  computation 
lengthy,  the  model  is  tractable  in  the  details.  The  greatest  difficulty 
is  in  accurately  calculating  the  effect  of  intervalley  scattering.  In 
this  analysis  a  method  due  to  Herring21  is  followed.  The  transition  of  an 
electron  from  one  valley  to  another  (in  k  space)  is  accompanied  by  a  large 
change  in  momentum.  The  phonons  which  participate  in  this  scattering 
process  lie  near  the  edge  of  the  Brillouin  zone.  At  the  Brillouin  zone 
boundary,  the  acoustic  and  optical  branches  are  either  degenerate  or  close 
to  one  another.  The  frequency  of  the  intervalley  mode  has  a  value 
between  the  frequency  of  the  two  modes  at  the  zone  boundary,  with  the 
relaxation  time  for  intervalley  scattering  of  carriers  given  by:- 

w2f[(EABei)+1]l/2]  +  exp(8i/T)Re[C(EABei)“i]^ 

Tf1=  L - 1 - 

[exp(0i/T)-l] 

where  W2  is  a  rate  constant.  Averaging  over  a  Maxwell-Boltzmann 
distribution  function  gives  for  thermal  carriers. 

<Ti~1>  =  w2 (2Z'/n)^  K1(Z')sinh(Z') 

where  is  a  modified  Bessel  function  of  the  second  kind  and  Z'  =  0^/2T. 

The  total  relaxation  time  of  carriers  is  obtained  by  adding  the 
inverse  relaxation  times  for  the  acoustic  scattering  Tac  and  intervalley 
scattering  r^.  The  ratio  nA0,  where  n  is  the  mobility  can  be 
calculated21  as  a  function  of  T/0^  and  the  mobility  obtained  for  various 
values  of  Wj/V^ ;  here  n0  =  Mac( T/0 i ) 3 ^ 2 ,  0^  =  ftw^AB  an^  is  the  rate 
constant  for  intervalley  scattering. 

The  use  of  Maxwell-Boltzmann  statistics  may  well  be  inappropriate  at 
the  high  carrier  concentrations  present  in  thermoelectric  semiconductors. 
However,  lead  telluride  exhibits  intervalley  scattering  to  a  relatively 
small  degree.  Consequently,  any  shortcomings  of  the  theoretical  model 
which  may  result  from  the  use  of  Maxwell-Boltzmann  statistics  will  have 
little  effect  on  the  calculated  values  of  the  electronic  transport 
coefficients.  The  calculation  of  the  Lorenz  factor  and  the  ratio  xe/)>L 
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are  based  upon  Ferai-Dirac  statistics  and  are  valid  at  high  carrier 
concentrations. 

4.  Thermoelectric  Figure  of  Merit 

The  dimensionless  thermoelectric  figure  of  merit  ZT  for  a  two  band 
model  can  be  written  as”:- 


(«y  -  c^on')1 

ZT  *=  - — 

(°n  +  °p  )(1+°n' rn +a'prp)  +  «n,®p'(Bn+Vtg>1 

where  a'  and  o'  are  the  dimensionless  forms  of  the  Seebeck  coefficient  and 
electrical  conductivity  . 

V  THEORETICAL  ANALYSIS 
1  Lead  Telluride 

(i)  Relative  reduction  in  thermal  conductivity  due  to  phonon-grain 
boundary  scattering. 

An  estimate  can  be  made  of  the  relative  reduction  in  the  lattice 
thermal  conductivity  of  unalloyed  lead  telluride  with  reduction  in  grain 
size  compared  with  single  crystal  material  using  the  model  outlined.  The 
results  of  this  calculation  are  displayed  in  Figure  1. 


Figure  1. 

Plot  xL(sintered)/Xr( single  crystal) 
for  unalloyed  lead  telluride  at  300K 
as  a  function  of  grain  size  and 
level  of  doping.  k0  *=  1.0,  A  =  0; 
curves:  1,  B  =  0.050:  2,  B  =  0.010; 
3,  B  «  0.005:  4,  B  *=  0. 


t  i^'l 
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Although  there  is  some  difficulty  in  making  an  accurate  prediction  of  the 
dependence  of  the  lattice  thermal  conductivity  on  the  various  parameters 
because  of  a  lack  of  experimental  data  and  the  uncertainty  in  the 
appropriate  value  of  k0  it  is  possible  to  estimate  the  range  over  which 
the  results  will  vary.  In  undoped  lead  telluride  the  reduction  in  lattice 


thermal  conductivity  is  given  by  { l-xLsintered/xLsingle  crystal)  x  100  for 
a  mean  grain  size  of  ljim  is  around  6  percent.  For  doped  lead  telluride 
(10*4m-#-10**m“*)  this  reduction  is  decreased  to  about  5  percent. 

(ii)  The  ratio  xe/xL  and  the  Lorenz  factor  t 


Of  central  importance  in  a  quantitative  evaluation  of  the  figure  of 
merit  is  the  contribution  to  the  thermal  conductivity  which  arises  from 
charge  carriers.  This  contribution  xe/x L  can  conveniently  be  expressed  in 
terms  of  the  quantity  o'*.  The  calculated  values  of  the  ratio  xe/XL  and 
the  Lorenz  factor  are  displayed  in  Figure  2  as  a  function  of  the  reduced 


Figure  2.  The  ratio  xe/xL  and  the 
Lorenz  factor  £  as  a  function  of 
reduced  Fermi  engergy  ({)  for  lead 
telluride  at  300K.  (a)  acoustic 

phonon  scattering  without 

intervalley  scattering  and  (b)  polar 
optical  scattering.  Curves  A  and  B 
refer  to  parabolic  and  non  parabolic 
bands  respectively. 


o  1 


Figure.  3.  The  ratio  xp/xL  as  a 
function  of  reduced  Fermi  energy  ({) 
for  lead  telluride  at  different 
temperatures.  A  non  parabolic 
energy  band  is  considered  with  Bg  = 
kgT/Eg;  acoustic  phonon  scatterng. 
Curves  (A)  300K  <B)  500K  (C)  700K 
(D)  900K. 


Fermi  energy  ({)  for  parabolic  and  non  parabolic  bands.  Scattering  of 
electrons  by  acoustic  and  optical  modes  have  been  included.  Important 
material  parameters  for  lead  telluride  are  collected  in  Table  1.  In 
Figure  3  is  displayed  the  ratio  Xg/Xj,  versus  {  for  lead  telluride  at 
different  temperatures  for  a  non  parabolic  energy  band.  The  temperature 
variation  of  Sg  is  obtained  from  the  relationship  =  k0T/Eg.  The  Lorenz 


Table  1.  Important  Material  Parameters  for  PbTe 


1.39  x  1011 
(Nm~2 } 


1.7 

(vm-1K-1 ) 


0.32  400 

(eV) 


factor  is  plotted  in  Figure  4  as  a  function  of  temperature  for  different 
carrier  concentrations. 


Figure  4.  The  Lorenz  factor  £  for 
lead  telluride  as  a  function  of 
temperature  at  different  carrier 
concentrations.  Acoustic  phonon 
scattering.  Curves  (A)  3xl02Sm"’, 
(B)  102Sm~\  (C)  5xl024m"3,  (D) 
1024m'5  (E)  2xl02 sm“5 . 


300  S00 


It  is  apparent  from  Figure  2  that  the  ratio  xe/xL  corresponding  to  a 
particular  value  of  <  is  considerably  reduced  when  the  effect  of  non 
parabolicity  is  taken  into  account.  In  thermoelectric  applications  the 
behavious  of  the  material  in  the  region  of  $opt  (value  of  the  reduced 
Fermi  energy  which  optimises  the  thermoelectric  figure  of  merit)  is  of 
considerable  interest.  In  Table  2  the  ratio  xe/xL  for  PbTe  corresponding 
to  $opt  at  different  temperatures  is  presented  for  parabolic  and  non 
parabolic  energy  bands.  Carrier  scattering  has  been  taken  to  be  by 
acoustic  phonons  only. 


Xe/XL 

(parabolic  band) 

(non-parabolic  band) 

Xe/XL* 

(non  parabolic  band) 

T(k) 

*opt 

6g=0 

8g=0.08 

Bg=lAg(T) 

300 

-0.75 

0.50 

0.27 

0.25  (0.15)* 

500 

-0.90 

0.70 

0.35 

0.30 

700 

-1.20 

1.00 

0.50 

0.40 

900 

f 

»-* 

o 

1.25 

0.65 

0.50 

r  obtained  from  calculation  of  the  figure  of  merit  and  corresponds 
to  about  2xl024m“s. 

r  (3  =0.08  is  a  realistic  estimate  of  0g  at  300K  but  because  is 
+  temperature  dependent  it  becomes  less  realistic  aty  higher 

temperatures .  ^^A*****^ 

?  Eg ( T )  =  Eg(  300 )  +  4  X  10-4  (T-300 ) ;  <g(T)  =@§7=  Eg(T)ABT. 
♦Includes  intervalley  scattering  W2Aj  =  0.5. 

Table  2.  The  ratio  Xe/Xr  for  lead  telluride  corresponding  to  at 

different  temperatures  and  appropriate  to  parabolic  and  non-parabolic 
bands;  acoustic  phonon  scattering. 

In  Table  3  are  collected  the  room  temperature  values  of  xe/xL  at 
optimum  doping  appropriate  to  optical  phonon  scattering  and  acoustic 
phonon  scattering,  separately  and  when  combined  (non-parabolic  energy 


bands ) . 

Xe/XL 

Xe/XL 

Xe/XL 

Xe/XL 

(acoustic) 

(optical) 

(combined) 

(observed) 

0.25(0.5)t 

0 . 50 ( 0 . 8 )  + 

0.17 

0.15 

t  obtained  using  parabolic  bands 

Table  J.  Room  temperature  values  of  the  ratio  xe/xL  at  optimum  doping 
(non-parabolic  energy  bands) 

(iii)  Electronic  Thermal  Conductivity 

The  inclusion  of  non-parabolicity  in  the  calculation  of  both 
scattering  mechanism  considerably  reduces  the  electronic  contribution  to 
the  thermal  conductivity.  The  calculated  value  of  the  combined  scattering 
is  in  very  good  agreement  with  the  experimentally  obtained  values23. 
Combining  the  effect  of  the  two  scattering  mechanisms  is  not 
straightforward.  A  rough  estimate  can  be  obtained  by  assuming  l  at  topt 
to  be  the  same  for  both  cases;  the  inverse  electrical  conductivities  can 


then  be  added  to  give  the  combined  values. 

In  Figure  5  are  displayed  the  calculated  values  of  the  electronic 


thermal  conductivity  of  n-type  lead  telluride 


Fipure  5.  Electronic  thermal 
conductivity  (xe)  for  n-PbTe  as  a 
function  of  carrier  concentration 
and  temperature  (acoustic 
scattering  '  and  non-parabolic 
bands).  Curves  1,2, 3, 4  correspond 
to  temperatures  of  300,  450,  600, 
750K,  respectively. 


figure  6.  Electronic  thermal  conductivity 

versus  carrier  concentration.  -  2x 
Curves  1,2, 3, 4  correspond  to  the 
temperatures  of  300,  450,  600,  750K,  E  20 

respectively.  *  , 


as  a  function  of  carrier  concentration  at  temperatures  of  300,  450,  600, 
and  750K  assuming  that  the  electrons  are  scattered  by  acoustic  phonons.  At 
carrier  concentrations  10Mm"!  the  electronic  contribution  becomes 
dominant;  consequently  the  contribution  due  to  polar  optical  scattering 
(Figure  6)  should  be  included  in  the  calculation  of  \e. 

(iv)  Total  Thermal  Conductivity 

A  reasonable  estimate  of  the  temperature  variation  of  the  total 
thermal  conductivity  can  be  obtained  by  assuming  a  temperature  dependence 
of  ‘/t  for  the  lattice  thermal  conductivity.  The  results  of  this 
calculation  are  plotted  in  Figure  7  assuming  acoustic  phonon  scattering. 
In  Figure  8  the  total  thermal  conductivity  is  plotted  as  a  function  of 
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temperature  for  two  different  carrier  concentrations;  n=5xl024cm"*  and 
102‘cm_s.  Curves  1  and  2  correspond  to  acoustic  phonon  scattering  and 
include  the  effect  of  non-parabolic  energy  bands.  Curves  1*  and  2'  are 
the  total  thermal  conductivity  with  xe  given  by  xe-1  =  xe-1  (acoustic)  + 
xe_1  (optical).  Evidently  the  inclusion  of  polar  optical  scattering  has 


Figure  7.  Total  thermal  conductivity 
x(=xL+\e)  as  a  function  of  carrier 


concentration  and  temperature  (acoustic 
scattering).  1,1'-300K;  2.2'-450K; 

3,3' -600K;  4,4'-750K.  (Dashed  lines: 

parabolic  bands;  solid  lines: 
non-parabolic  bands. 


Figure  S.  Total  thermal  conductivity 
(x)  versus  temperature.  1,1 '-5  x 
I024m“3;  2,2' -102  5m”J ;  l',2' 

include  scattering  of  electrons  by 
polar  optical  modes  along  with  the 
acoustic  modes.  Experimental 
points  by  Efimova  et  al,  [13]: 

A5xl02  4m“3  and  ■102  5nfs 


3°  y.A 

?  2.6  *\\ 
F  *.\  > 


^  2.2  h  r\\  ' 


300  tOO  500  600  700 
T  Ik) 

an  appreciable  effect  on  the  total  thermal  conductivity  at  carrier 
concentrations  in  excess  of  5xl024m"3,  The  calculation  of  the  lattice 
thermal  conductivity  in  Figure  8  also  takes  into  account  the  scattering  of 
phonons  by  free  electrons.  The  parameter  B  was  adjusted  to  the  values  of 
0.01  and  0.02,  curves  (1,1')  and  (2,2')  respectively.  Excellent  agreement 
is  obtained  between  the  calculated  values  and  reported  experimental  data 
at  around  room  temperature24.  However,  it  is  apparent  that  the  slopes  of 
the  theoretical  curves  are  somewhat  greater  than  those  based  upon 
experimental  data.  The  inclusion  of  a  photon  contribution  to  the  thermal 
conductivity  at  higher  temperatures  is  likely  to  be  small  for  doped 
samples19  and  is  unlikely  to  account  for  the  differences  in  slope.  In 
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Figure  9  is  displayed  the  effect  of  phonon-grain  boundary  scattering  on 
the  room  temperature  thermal  conductivity  of  small  grain  size  lead 
telluride.  Earlier  calculations  of  this  effect* *  did  not  take  into 
account  the  electronic  contribution  to  the  thermal  conductivity. 


Figure  9.  Total  thermal  conductivity 
versus  carrier  concentration  for 
"single  crystal"  lead  telluride  and 
fine-grained  material  with  mean 
grain  size  lfim. 


The  calculation  assumes  an  absence  of  disorder  (A^o)  and  consequently  the 
reduction  in  thermal  conductivity  corresponds  to  the  situation  when 
boundary  scattering  is  relatively  less  effective  in  scattering 
phonons. 17'*‘  In  undoped  lead  telluride  with  a  mean  grain  size  of  l»m  the 
reduction  in  thermal  conductivity  will  be  about  6%  whereas  in  the  range  of 
optimum  doping  encountered  in  thermoelectric  applications  (n"2xlO*4m“*  in 
lead  telluride)  this  reduces  to  about  5%. 

2.  Disordered  Lead  Telluride 
(i)  Introduction 

As  mentioned  previously,  phonon-grain  boundary  scattering  will  be 
enhanced  in  alloys.  The  initial  theoretical  model,  which  was  developed 
to  obtain  an  estimate  of  the  relative  reduction  in  the  lattice  thermal 
conductivity  of  lead  telluride  with  decrease  in  grain  size,  was  extended 
to  include  disordered  lead  telluride  (alloys).  An  estimate  was  obtained 
of  the  dependence  of  the  thermal  conductivity  and  thermoelelectric  figure 
of  merit  on  grain  size,  level  of  doping  for  a  highly  disordered  alloy. 
Preliminary  calculations  indicated  that  calculations  involving  relative 
changes  in  the  thermoelectric  figure  of  merit  are  fairly  insensitive  to  the 
inclusion  of  refinements  in  the  theoretical  model.  A  two  band  model  with 
a  parabolic  multivalued  structure  was  considered,  acoustic  scattering  was 
taken  as  the  dominant  scattering  mechanism,  intervalley  scattering  was 
neglected  and  no  distinction  made  between  conductivity  effective  mass  and 
the  density  of  states  effective  mass. 
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(ii)  Reduction  in  the  lattice  thermal  conductivity 

Plots  of  the  ratio  xsintered/xsingle  crystal  at  room  temperature  for 
highly  disordered  alloys  of  lead  telluride  (A=5)  is  displayed  in  Figure  10 
as  a  function  of  grain  size  and  level  of  doping;  in  material  with  a  mean 
grain  size  of  about  l*m  the  reduction  in  lattice  thermal  conductivity 
compared  to  equivalent  single  crystal  material  is  in  the  range  11-13 
percent.  ,  (um) 


Figure  IQ.  Plot  of 
xL( sintered) xL( single 
crystal  )fi  |#  for  a  highly 
disordered  alloy  of  lead 
telluride  of  grain  size  and 
level  of  doping.  kg 
=1.0.A=5.0;  curves 
1,B=0.050;  2,B=0.010;  3, 

B=0.005;  4,  B=0 
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Figure  11.  Plots  of  the 
thermoelectric  figure  of  merit  Z 
for  a  highly  disordered  alloy  of 
lead  telluride  with  three  different 
carrier  concentrations  n,  as  a 
function  of  temperature:  A, 
n  =  5xl02  4m~3 ;  B,  n  =  1025ra~3; 

C,  n  =  2xl025nf3.  Curves:  - xL 

corresponds  to  material  with  a  mean 
grain  size  -Imi®;-  -  -  -xL 
corresponds  to  'single  crystal'  or 
large  grain  size  material 
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( i i i  )  Effect  of  small  grain  size  on  the  thermoelectric  figure  of  merit . 

As  indicated  in  Table  2  preliminary  calculations  of  the 
thermoelectric  figure  of  merit  (Z)  of  unalloyed  lead  telluride  indicated 
that  Z  was  optimised  at  room  temperature  at  a  reduced  Fermi  Energy  ( { )  of 
-0.75;  which  corresponded  to  a  carrier  concentration  of  about  2xl024m"3. 
The  thermoelectric  figure  of  merit  for  a  highly  disordered  alloy  of  lead 
telluride  at  their  different  carrier  concentrations  around  optimum  doping 
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is  shown  in  Figure  11  as  a  function  of  temperature.  The  temperature 
dependence  of  the  effective  mass  m*  and  the  energy  gap  Eg  is  taken  into 
account  using  the  relationship 

1  dm*_  1  dEg 
m*  dT  ^  ~d!I 

dE 

Eg(T)  =  Eg{T0)  -  (T-300) 


with  T0  =  300K,  dEg/dT  =  4xl0“4,  m*  =  0.2mo  at  300K,  XL  at  300K  is 
l.TOwnT1^1  and  is  assumed  to  vary  inversely  with  temperature.  The 
variation  of  carrier  concentration  n  with  temperature  is  obtained  by 
relating  the  variation  of  Seebeck  coefficient  with  reduced  Fermi 
potential26  to  the  temperature  dependence  of  the  Seebeck  coefficient  for 
various  carrier  concentrations. 2*  In  Figure  11  a  comparison  is  drawn 
between  "single  crystal"  and  small  grain  size  material.  It  is  evident 
that  the  thermoelectric  figure  of  merit  of  highly  disordered  lead 
telluride  type  material  with  a  mean  grain  size  of  -l^m  is  about  10  percent 
higher  than  equivalent  "single  crystal"  or  large  grain  size  material. 

3 .  Alloys  based  upon  lead  telluride 
(i)  Introduction 

A  realistic  theoretical  model  for  lead  telluride,  which  included  a 
multivalued  structure,  intervalley  and  intravalley  scattering,  non 
parabolic  energy  bands,  acoustic  phonon  and  optical  phonon  scattering;  has 
been  developed  and  close  agreement  obtained  between  theoretical  and 
measured  transport  properties.  This  model  was  then  employed  in 
identifying  alloys  based  upon  lead  telluride  with  the  highest  potential 
for  improvement  in  their  figure  of  merit  as  a  result  of  phonon-grain 
boundary  scattering.  PbSnTe  and  PbGeTe  are  particularly  favoured  because 
large  differences  in  atomic  masses  of  the  constituent  atoms  give  rise  to 
substantial  alloy  disorder  scattering  and  these  materials  were 
investigated. 

{ i i )  Reduction  in  the  lattice  thermal  conductivity  of  PbSnTe  and  PbGeTe. 

The  results  of  calculating  the  ratio  xL(sintered)/xL(single  crystal) 
as  a  function  of  grain  size  L,  conveniently  expressed  in  terms  of  a 
parameter  D  where  D  =  CT  and  T  is  the  absolute  temperature  for  PbSnTe  and 
PbGeTe  are  displayed  in  figures  12  and  13. 
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Figure  12.  x,  the  ratio 
Vr ( sintered )/xL( single  crystal), 
plotted  as  a  function  of  the 
parameter  D  and  grain  size  L  for 
PbTe-SnTe  at  300K.  Curve  I,  B=0 
(undoped);  curve  II,  B=0.01 
(optimally  doped). 


Figure  13.  x,  the  ratio 
Xt (sintered)/xL( single  crystal), 
plotted  as  a  function  of  the 
parameter  D  and  grain  size  L  for 
PbTe-GeTe  at  300K.  Curve  I,  B=0 
(undoped);  curve  II,  B=0.01 
(optimally  doped). 
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The  basic  characteristics  of  the  PbTe  system  has  been  retained  in  the 
calculation  and  the  effect  of  alloying  with  SnTe  and  GeTe  is  to  produce 
disorder  which  can  effectively  scatter  the  high  frequency  phonons.  This 
has  the  effect  of  enhancing  the  effectiveness  of  phonon-grain  boundary 
scattering.  Although  there  is  some  difficulty  in  making  accurate 
predictions  of  the  dependence  of  the  lattice  thermal  conductivity  on  the 
various  parameters  because  of  a  lack  of  experimental  data  and  an 
uncertainty  in  fixing  an  appropriate  value  of  k0  it  is  possible  to 
estimate  the  range  o-er  which  the  results  will  vary. 

It  is  concluded  that  in  undoped  PbSnTe  and  PbGeTe  (B=0)  with  a  mean 
grain  size  of  0.5^  at  300K  the  percentage  reductions  in  lattice  thermal 
conductivity  given  by  [1  -  xL{sintered)/xL((single  crystal))  x  100  are  15 
percent  and  18  percent  respectively.  In  optimally  doped  material 
(assuming  carrier  concentration  to  be  similar  to  that  of  unalloyed  lead 
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telluride  the  reductions  are  11  percent  and  14  percent  respectively.  A 
further  reduction  in  mean  grain  size  to  0.25/xm  would  decrease  the  lattice 
thermal  conductivity  of  lead  telluride  by  17  percent  and  PbGeTe  by  21 
percent.  This  is  approaching  the  limit  of  the  beneficial  effect  of  a 
reduction  in  grain  size  as  in  this  region  electron  grain  boundary 
scattering  becomes  significant  and  leads  to  an  undesirable  decrease  in  the 
electrical  conductivity. 

VI  EXPERIMENTAL  PROGRAMME  OF  WORK. 

1.  Introduction 

The  theoretical  model  outlined  in  Section  I  provided  a  guide  on  the 
reduction  in  lattice  thermal  conductivity  with  decrease  in  grain  size.  The 
objective  of  the  experimental  programme  of  work  was  to  substantiate  the 
theoretical  predictions  and  towards  this  end  it  was  decided  *o  attempt  to 
prepare  compacted  material  with  a  mean  grain  size  of  <5jzm  and  <0.5pim. 

2.  Charge  material  preparation 


(i)  <5um  grain  size 

The  starting  material  used  in  the  investigation  was  a  pulled,  large 
grain  size  (up  to  0.5cms)  ingot  of  2N-lead  telluride  supplied  by  Global 
Thermoelectrics*  2N  is  a  3M*  designation  for  lead  telluride  to  which 

approximately  0.03  molelcular  percent  Pbl2  has  been  added  to  increase  the 
carrier  concentration.  The  ingot  was  crushed  under  methanol  in  an  agate 
pestle  then  ground  wet  for  one  hour  in  a  two  ball  vibromill.  The  powder 
was  sieved  through  a  5jim  British  Standard  microsieve  using  methanol  as  a 
vehicle  and  assisted  by  ultrasonic  vibrations. 

(ii)  <0.5um  grain  size 

Decreasing  the  grain  size  of  the  charge  material  by  an  order  of 
magnitude  prescribed  handling  problems.  Submicron  size  particles  are 
airborne  at  low  velocity  draughts  and  all  handling  stages  of  the  powder 
were  carried  out  when  wet  with  methanol.  The  <5*zm  grain  size  sieved 
fraction  obtained  as  described  above  was  ground  in  an  agate  ball  mill  for 
one  hour  using  methanol  as  a  vehicle.  A  photomicrograph  of  the  resultant 
powder  dispersed  in  NONIDET  P42  and  viewed  in  transmitted  light  is  shown 
in  Figure  14.  The  powder  rests  on  the  base  of  a  5^m  aperature  microplate 

*  Global  Thermoelectric  Power  Systems  Ltd 
PO  Box  400 
Bassano 

Alberta,  Canada 
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Figure  14.  Photomicrograph  of 
<0.5/zm  lead  telluride  powder 
dispersed  in  NONIDET  P42;  5nm 
microsieve  operature  shown  as  a 
size  comparison. 
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sieve  and  the  operative  provides  a  convenient  size  comparison  for  the 
particles.  Considerable  difficulty  was  encountered  in  satisfactorily 
dispersing  these  very  small  particles  even  when  assisted  by  ultrasonic 
vibrations  as  is  evident  from  the  agglomerated  appearance.  No  attempt  was 
made  to  measure  the  size  of  the  particles,  but  evidently  the  vast  majority 
of  particles  are  less  than  0.5/im.  The  charge  is  introduced  into  the  die 
as  a  "slurrey"  and  the  methanol  evaporated  off  insitu. 

3.  Hot  Press 

It  is  reported  in  the  literature2 7  that  lead  telluride  can  be 
successfully  compacted  using  hot  or  cold  pressing  techniques.  However, 
lead  telluride  is  known  to  possess  a  relatively  high  vapour  pressure  and 
losses  of  constituents  result  in  a  change  in  transport  properties; 
consequently  initial  attempts  were  made  to  cold  compact  the  powders.  The 
arrangement  employed  in  this  present  study  is  shown  in  Figure  15.  A  cold 
compaction  pressure  of  300  MPa  is  recommended,  consequently  emphasis  was 
placed  on  strength  and  stability  with  the  die  and  supporting  column 
machined  in  a  number  of  interlocking  sections.  Provision  for  rapidly 
heating  the  die  is  made  by  locating  the  die  inside  a  relatively  large 
graphite  cylinder.  The  die  and  plungers  are  frabricated  from  TZM  alloy 
and  machined  to  a  press  fit.  The  inside  surface  of  the  die  is  painted 
with  liquid  graphite  and  the  faces  of  the  plungers  separated  by  high 
density  graphite  spacers. 
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Hydraulic  Pressure 

Figure  15. 

Hot  Press  Arrangement 


(a)  Pt/Pt-Rd  shielded  thermocouple  (b)  Silver  steel  packing  block 

(c)  Graphite  susceptor  (d)  TZM  die 

(e)  Mild  steel  support  (f)  silica  wool  thermal  insulation 

(g)  r.f  induction  coil  (h)  PTFE  insulator 

(i)  Water  cooling 

4 .  Pressing  procedure 

A  series  of  cold  compactions  were  carried  out  on  <5jzm  powder  over  a 
wide  range  of  pressures.  Single  crystal  density  is  8.25  gm  cm"3  and  the 
density  of  the  compacts  increased  with  pressing  pressure  from  5.0  gm  cm"3 
at  150  MPa  to  7.80  gm  cm"3  at  400  MPa  as  shown  in  Figure  16.  Although  it 
should  be  noted  that  this  graph  is  only  a  guide  as  its  shape  is  strongly 
dependent  upon  the  state  of  the  particle  surface  and  particle  size. 


Figure  16.  Dependence  of  density 
of  PbTe  compacts  with  pressing 
pressure . 


Cracks  develope  in  compacts  pressed  at  >300MPa  and  they  invariably  break 
into  sheet-like  pieces  in  planes  perpendicular  to  the  pressing  direction. 
The  application  of  heat  greatly  reduces  the  requires  compaction  pressure. 
Mechanically  strong  compacts  with  a  grain  size  of  <o.5/zm  and  densities  of 
around  8.1  gm  cm-3  were  prepared  by  hot  pressure  sintering  for  10  minutes 
at  a  pressure  of  approximately  100  MPa  and  a  temperature  of  HOOK. 
Compacts  of  similar  densities  were  also  prepared  with  a  grain  size  of 
<0.5jm,  although  these  compacts  were  not  as  mechanically  strong.  As  the 
main  objective  of  this  stage  of  the  research  programme  was  to  determine 
the  change  in  thermal  properties  with  grain  size  rather  than  evaluate  the 
effect  of  high  temperature  compaction  on  the  electrical  properties  of  PbTe 
subsequent  physical  characterisation  and  measurement  of  transport 
properties  were  made  on  hot  pressed,  high  density  compacts. 


Figure  17.  Photomicrograph  of 
<0.5/*m  grain  size  compacted 
lead  telluride. 


5.  Physical  Properties 

The  surface  of  the  high  density,  pressed  discs,  after  polishing  down 
to  ‘/up  size  is  smooth  and  void  free  when  examined  by  optical 
microscopy.  The  grain  structure  of  the  compact  is  revealed  by  heating  in 
an  iodine  etch  at  368K  for  5  minutes  (10H20,  5gm  NaoH  and  0.2gm  I2)28.  A 
photomicrograph  of  <0.5*iro  grain  size  compacted  material  is  shown  in  Figure 
14.  It  is  apparent  that  the  "larger  grains"  are  agglomerates  and  that 
little,  if  any,  grain  growth  has  taken  place.  Density  of  the  compacts 
were  determined  by  the  method  of  hydrostatic  weighing.  Error  in  density 
determination  is  less  than  0.5  percent. 

6.  Transport  Properties. 

(i)  Introduction 

The  transport  properties,  and  in  particular  the  electrical  properties 
are  determined  by  the  stoichiometric  ratio  of  its  components.  An  excess 
of  lead  induces  n-type  electrical  behaviour  and  an  excess  of  tellurium, 
p-type  by  behaviour,  carrier  concentration  of  about  3xl023m“3  can  be 
induced  in  this  way.  In  thermoelectic  application  carrier  densities  an 
order  of  magnitude  larger  are  required  than  those  which  can  be  achieved  by 
changing  the  stoichiometric  ratio.  These  higher  values  of  carrier 
concentration  are  obtained  by  introducing  foreign  molecular  species  such 
as  Pbl2  (elecrons)  and  Na  (holes).  The  solid  solubility  limit  of  the 
components  are  temperature  dependent  and  is  determined  by  the  annealing 
temperature.  The  carrier  concentration  also  depends  upon  the  rate  of 
quenching  with  more  carriers  remaining  in  solution  if  the  material  is 
rapidly  quenched. 

(ii)  Seebeck  coefficient  and  electrical  resistivity  measurements. 


Seebeck  coefficient  measurements  were  made  using  a  hot  probe? 
accuracy  of  measurement  is  ±  3  percent.  Electrical  resistivity 
measurements  were  made  using  the  four  probe  method;  accuracy  ±  2  percent. 
Measurements  of  Seebeck  coefficient  and  electrical  resistivity  were  made 
on  the  pressed  compacts  and  compared  with  values  obtained  from 
measurements  on  discs  having  identical  geometry  and  cut  from  the  original 
"single  crystal"  ingot.  The  Seebeck  coefficient  and  electrical 
conductivity  of  "as  compacted"  material  differs  considerably  from  that  of 
single  crystal.  High  temperature  annealing  is  an  established  procedure 
for  stabilizing  the  transport  properties  of  semiconductors  and  the 
compacted  materials  have  been  subjected  to  a  variety  of  annealing 
sequences.  The  effect  of  rapidly  cooling  the  sample  (quenching  in  water) 
is  to  "freeze  in"  the  large  number  of  carriers  present  at  the  higher 
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temperature  and  hence  reduce  the  Seebeck  coefficient.  Rapid  cooling, 
however,  is  accompanied  by  an  increase  in  the  strain  energy  of  the  sample, 
and  results  in  a  decrease  in  the  carrier  moblity  and  consequent  increase 
in  resistivity.  To  date  electrical  properties  closest  to  single  crystal 
have  been  obtained  by  annealing  the  compact  in  an  argon  atmosphere  at 
1040K  for  3  hours  followed  by  4  hours  at  700K  and  slowly  cooling  to  room 
temperature.  Seebeck  coefficients  very  close  to  that  of  single  crystal 
have  been  obtained  whilst  the  electrical  conductivity,  although 
substantially  reduced,  remains  significantly  higher  than  single  crystal 
value. 

(iii)  Thermal  diffusivity  and  thermal  conductivity  measurements. 

Room  temperature  thermal  diffusivity  measurements  were  made  on  hot 
pressed  samples  using  a  laser  flash  technique29.  Discs  5mm  in  diameter 
were  cut  from  the  <5jm  compacted  material  using  an  ultrasonic  drill. 
Thermocouple  contacts  to  the  rear  of  the  sample  are  made  using  pressure 
contacts.  Optimum  sample  thickness  for  thermal  diffusivity  measurements 
was  about  lmm,  but  samples  of  this  thickness  were  too  fragile  and 
invariably  fractured  when  subjected  to  the  contacting  pressure  which  was 
necessary  to  provide  low  electrical  resistance  contacts.  Consequently 
samples  2mm  thick  were  employed.  As  indicated  previously  compacts  of 
<0.5/xra  grain  size  material  were  relatively  weak  and  proved  difficult  to 
machine.  Thermal  conductivity  measurements  were  made  on  the  polished 


Figure  IS.  The  reduction  in 
thermal  conductivity  of  PbTe 
with  decrease  in  grain  size 

(L)  -  optimally  doped . 

undoped  (theoretical  curves) 

•  Experimental  values. 
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surface  of  these  very  small  grain  size  samples.  A  suitable  disc 
thickness  is  lcra.  A  thermal  comparator  method  was  employed  and  again  a 
comparison  made  between  the  compacted  sample  and  one  of  identical  geometry 


cut  from  the  "single  crystal"  ingot.  The  reduction  in  thermal 
conductivity  with  decrease  in  grain  size  (at  room  temperature)  is  shown  in 
Figure  15. 

vii  Discussion  and  conlusion 

All  principle  objectives  of  the  programme  of  research  have  been 
achieved.  A  semiquantitative  theoretical  model  of  lead  telluride  has  been 
used  to  estimate  the  relative  reduction  in  lattice  thermal  conductivity, 
compared  to  that  of  single  crystal,  which  accompanies  the  use  of  small 
grain  size  material.  It  was  predicted  that  in  material  with  a  mean  grain 
size  of  14m,  the  reduction  would  be  around  5  percent.  A  procedure  has  been 
developed  for  comminuting  very  small  grain  size  material  and  a  number  of 
high  density  compacts  of  lead  telluride  successfully  prepared. 
Measurements  on  small  grain  size  compacts  substantiated  the  predicted 
reduction  in  thermal  conductivity  with  decrease  in  grain  size. 

A  realistic  theoretical  model  has  been  developed  for  lead  telluride. 
Very  good  agreement  was  obtained  between  the  theoretical  values  of  the 
thermal  conductivity  and  experimental  values  cited  in  the  literature.  The 
thermoelectric  figure  of  merit  is  optimised  at  a  carrier  concentration  of 
about  2xlOz  4m“3 .  Provided  the  electrical  properties  can  be  maintained 
close  to  "single  crystal"  value  the  thermoelectric  figure  of  merit  in 
highly  disordered  alloys  of  lead  telluride  with  a  mean  grain  size  of  -ljim 
and  optimum  doping  is  calculated  to  be  about  10  percent  higher  than 
equivalent  single  crystal  values. 

The  model  has  been  used  to  explore  the  thermal  properties  of  alloys 
based  upon  lead  telluride.  PbSnTe  and  PbGeTe  hold  out  the  best  potential 
for  improvement  in  the  figure  of  merit,  as  a  result  of  a  decrease  in  the 
lattice  thermal  conductivity  due  to  phonon-grain  boundary  scattering.  At 
room  temperature  the  lattice  thermal  conductivity  of  PbSnTe  and  PbGeTe, 
optimally  doped  and  with  a  grain  size  of  0.254m,  would  be  reduced  by  17 
percent  and  21  percent  respectively. 

Although  not  strictly  part  of  this  current  research  contract,  the 
behaviour  of  the  electrical  transport  properties  of  fine  grained  compacts 
is  very  relevant  to  any  attempt  to  improve  the  thermoelectric  figure  of 
merit.  The  results  of  a  limited  programme  of  work  in  this  area  indicate 
that  Seebeck  coefficients  very  close  to  single  crystal  values  can  be 
obtained  by  annealing  the  fine  grain  compacted  material.  Electrical 
resistivity  values  can  also  be  substantially  reduced  although  values 
remain  significantly  higher  than  equivalent  single  crystal  values.  A 
systematic  investigation  is  required  to  establish  the  annealing  procedure 


which  optimises  the  electrical  properties  of  these  compacts. 

It  is  concluded  that  the  thermal  conductivity  of  lead  telluride  type 
material  can  be  substantilaly  reduced  and  hence  the  thermoelectric  figure 
of  merit  significantly  improved  through  the  use  of  very  small  grain  size 
material,  provided  the  electrical  properties  can  be  maintained  at  values 
close  to  those  of  equivalent  single  crystal  material.  As  annealing 
procedures  are  routinely  employed  in  controlling  the  electrical  properties 
of  established  thermoelectric  materials  it  is  not  anticipated  that  the 
development  of  a  suitable  annealing  procedure  for  very  fine  grained  lead 
telluride  type  material  would  present  great  difficulties. 
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by  D*CT  where  T  is  the  temperature.  AB0  corresponds  to  unalloyed  material 
with  no  disorder  present,  B=0  corresponds  to  undoped  material  and  C=0  to 
single  crystal  material.  In  general  x(A,B,C=0)  *  Single  represents  the 
lattice  thermal  conductivity  of  a  doped  single  crystal  alloy  while 
x  (A,B,C)  *  Wintered  represents  the  lattice  thermal  conductivity  of  a 
compacted  (sintered)  alloy. 

3.  Electronic  contribution  to  the  thermal  conductivit 


(i)  Introduction 

Lead  telluride  and  its  alloys  possess  a  narrow  energy  band  gap  (Eg). 
These  materials  have  small  electron  effective  masses  and  the  density  of 
states  which  varies  as  m*H  is  also  relatively  small.  Consequently  a 
relatively  small  number  of  carriers  will  be  sufficient  to  fill  the  bands 
up  to  high  energy  levels.  The  effective  mass  of  electrons  near  the  top  of 
the  conduction  band  become  energy  dependent1 8 •  In  the  model  adopted  a  two 
band  conduction  band  is  considered  with  the  following  main  features;  the 
band  extrema  for  the  conduction  and  valence  bands  are  assumed  to  be 
located  at  the  same  k  value.  The  energy  separation  from  the  other  bands 
at  this  k  value  is  greater  than  the  main  energy  gap  and  the  momentum 
operator  has  non-zero  matrix  elements  between  the  states  corresponding  to 
the  extremal  points.  Non-parabolicity  and  the  effective  mass  of  both 
bands  are  defined  by  the  interaction  of  the  electron  and  hole  bands  only. 
The  energy  dispersion  law  is  of  the  form:- 


fi’kj2 


ft’k.f 


*  E  (1  +|  > 


k,  and  kn  are  the  transvei'se  and  longitudinal  components  of  the  wave 
vector  and  m*0  and  m*0  ^np  the  components  of  the  effective  mass  tensor 
near  the  band  extremum.  This  is  usually  referred  to  as  the  "Kane  model". 
The  energy  dependence  of  the  effective  mass  is  given  by 
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where  L^A^C) 
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here  x  -  hu/kBT  and  k0  «=  tu/tn  is  the  ratio  of  the  phonon  relaxation  tines 
for  Umklapp  and  Normal  processes,  kB  is  Boltzmann's  constant  and  T  is  the 
temperature.  There  is  no  advantage  in  this  present  study  to  undertake 
detailed  calculations  of  the  contributions  to  the  thermal  conductivity 
from  separate  polarizations  and  an  average  acoustic  branch  is  assumed. 
The  various  phonon  relaxation  times  are  as  described  in  the 
literature;14  18  17  the  phonon-electron  scattering  relaxation  time  is 
given  by 


B'x 


A,  B  and  C  are  given  by 
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Alloy  disorder 
scattering 


B  = 


2n! v0ab' 
kBi*V(i+ 


Free  carrier 
scattering 


The  various  parameters  have  the  same  meaning  as  in  the 
literature14  16  17  •  It  is  usual  to  express  C  in  terms  of  a  parameter  D, 
which  is  inversely  proportional  to  the  grain  size  L  and  they  are  related 


The  appropriate  expression  for  T(n)  which  includes  the  energy  dependence 
of  M  is  given  by**'10 

8Bgn(14egr>)' 

3(1+2 Bgn)* 

where  Cj  is  the  elastic  constant  related  to  the  average  sound  velocity, 
(p/Nv)  is  the  density  of  state  in  a  given  valley  and  c  is  the  deformation 
potential  constant. 

The  Lorenz  factor  is  given  by:- 

1  =  (’U./Vj-  6’ 
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where  6  =  1Lll/0LL2,  and  nLj  are  generalised. 

Fermi  integrals  and  defined  as:- 

nL|  (U  =  [-  nn[n(l+Bgn)]In(l+2Bgn)id^J 

The  reduced  electrical  conductivity  is  given  by:- 
c'  =  khu°Li 


V  2  - 1 - 

( *  L ) 

Kith  the  carrier  concentration  n  related  to  the  reduced  Fermi  energy 
by  the  expression:- 
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The  carrier  mobility  vac  due  to  acoustic  phonon  scattering  is  given 


by:- 
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3 .  ( i i i  )  Carrier  s catte ri ng  by  polar  optical  phonons 

The  carrier  relaxation  time  for  scattering  by  polar  optical  phonons 
is  given  by:0 
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The  Lorenz  factor  is  given  by 


t  =  [JLi2/°Lij)-6J 
where  6  =  iLli/°Lli 

In  the  calculations,  screening  effects  have  not  been  taken  into  account. 
Polar  optical  scattering  although  appreciable  in  lead  telluride  materials, 
is  not  the  dominant  scattering  mechanism.  Consequently  neglecting 
screening  effects  will  not  substantially  affect  the  overall  electronic 
transport  properties. 

The  reduced  electrical  conductivity  can  be  expressed  as:- 
o’  -  NyTJf(e0/T)md*°Li2/(XLmc*) 


with  K* 


where  f0  and  <„  are  the  static  dielectric  constant  and  the  dielectric 
constant  at  high  frequencies  respectively.  60  is  the  temperature  which 
corresponds  to  the  optical  phonon  frequency  at  the  zone  boundary. 

The  carrier  mobility  *z0p  due  to  optical  phonon  scattering  is  given 


by:- 
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where  K 
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Xe/XL 

(parabolic  band) 

v.1 

(non-parabolic  band) 

(non  parabolic  band) 

T(k) 

*opt 

Bg'0 

Bg-0.08 

3gBl/(g(T) 

300 

-0.75 

0.50 

0.27 

0.25  (0.15)* 

500 

-0.90 

0.70 

0.35 

0.30 

700 

-1.20 

1.00 

0.50 

0.40 

900 

-1.30 

1.25 

0.65 

0.50 

t  obtained  from  calculation  of  the  figure  of  merit  and  corresponds 
to  about  2xl0J4m”s. 

t  Bge0.08  is  a  realistic  estimate  of  6„  at  300K  but  because  B_  is 
+  temperature  dependent  it  becomes  less  realistic  atg  higher 

temperatures.  ttotUL' 

9  Eg (T )  *  Eg( 300 )  +  4  X  10"4(T-300);  <g(T)  *  Eg(T)/kBT. 

‘Includes  intervalley  scattering  V2/V)i  *  0.5. 

Table  2.  The  ratio  xe/Xt  for  lead  telluride  corresponding  to  (0Dt  at 
different  temperatures  ana  appropriate  to  parabolic  and  non-partfoblic 
bands;  acoustic  phonon  scattering. 

In  Table  3  are  collected  the  room  temperature  values  of  xe/xL  at 
optimum  doping  appropriate  to  optical  phonon  scattering  and  acoustic 


phonon  scattering, 
bands) . 

separately  and 

when  combined 

(non-parabolic  energy 

Xe/XL 

Xe/XL 

Xe/XL 

Xe/XL 

(acoustic) 

(optical) 

(combined) 

(observed) 

0.25(0.5)T 

0 . 50 ( 0 . 8 )  + 

0.17 

0.15 

t  obtained  using  parabolic  bands 

Tabic  3.  Room  temperature  values  of  the  ratio  xe/xL  at  optimum,  doping 
(non-parabolic  energy  bands) 

(iii)  Electronic  Thermal  Conductivity 

The  inclusion  of  non-parabolicity  in  the  calculation  of  both 
scattering  mechanism  considerably  reduces  the  electronic  contribution  to 
the  thermal  conductivity.  The  calculated  value  of  the  combined  scattering 
is  in  very  good  agreement  with  the  experimentally  obtained  values1 3 . 
Combining  the  effect  of  the  two  scattering  mechanisms  is  not 
straightforward.  A  rough  estimate  can  be  obtained  by  assuming  t  at  t0pt 
to  be  the  same  for  both  cases;  the  inverse  electrical  conductivities  can 


